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We analyze statistical properties of the separate multipole moments of the CMB temperature 
maps and find that the distribution tails are slightly non-Gaussian. Moreover, the deviation from 
Gaussianity peaks sharply at around / ~ 45 ± 10. If the detected non-Gaussianities should be 
attributed to the remaining foreground contamination from the galactic plane, an appearance of a 
similar peak in different frequency bands is unexpected. The presence of the peak was also confirmed 
by the analysis of weighted multipoles which is less susceptible to the galactic noise. A separate 
analysis of the multipoles with odd and even I + m shows two opposite types of non-Gaussianities 
which cannot be explained with only galactic foregrounds. We argue that cosmic strings or correlated 
gravitational waves could lead to the observed effect. 



I. INTRODUCTION 

The CMB radiation is known to be Gaussian to a very high degree, which is in a good agreement with predictions of 
the simplest model of a single scalar field slow- roll infiation llUa| . It is also known that considerable non-Gaussianities 
could develop in the models of multi-field infiation 0, 0, Such scenarios usually predict non- vanishing three- 

point correlation functions which are often discussed in the literature in the context of WMAP data analysis 0, [3, 
BE, [HI- At the same time the bispectrum is almost completely blind to non-Gaussianities from more complicated 
non-linear sources such as cosmic strings. Indeed, if the probability distribution of different multipole moments is 
nearly symmetric, then three-point correlation functions would not be very useful even if the distribution tails deviate 
significantly from a Gaussian. To detect such non-Gaussianities it seems necessary to use higher-point correlation 
functions (e.g. trispectrum). 

In this article, we discuss some methods which could be used to study the Gaussianity of the distribution tails. The 
analysis of WMAP temperature maps was performed for different multipole moments in the range 20 < Z < 100, where 
the data are cleanest but the number of realizations of each random variable is sufficient to study the probability 
distributions of different multipoles without having to worry too much about cosmic variance. Moreover, the separate 
aim's are (nearly) independent random variables which simplifies the analysis significantly. 

The paper is organized as follows. In Section II we describe the tests of Gaussianity which are employed to study the 
WMAP temperature maps. In Section III we discuss two different mechanisms involving cosmic strings and correlated 
gravitational waves, which could lead to the detected non-Gaussianities. 

II. DATA ANALYSIS 

It is often assumed that separate multipole moments of the CMB sky are Gaussian random variables. If so, then 
the WMAP temperature maps are completely specified by only a single function, the angular power spectrum 



C/ = (|a,„/). (1) 



Indeed, this would be true if the CMB temperature fluctuations were produced by one or few uncorrelated Gaussian 
processes. Then for a flxed (and large enough) I the values of ajm's must be described fairly well by a two dimensional 
Gaussian distribution with variance ax = = Generally, we expect all of I random complex variables^ on 



* Electronic address: | vitaly. vanchurin@physik.uni-muenchen.de| 

^ For the temperature fluctuations to remain real the condition = a*_„ must be satisfied. This means that strictly speaking we are 
dealing with 2Z -I- 1 random real numbers. One could get twice as much information by testing Gaussianity of the components of complex 
amplitudes separately, but to avoid unnecessary complications we restrict our analysis to multipole moments with m > 0. 
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Figure 1: l{l + l)G/(27r) (solid blue) and l{l + l)C/^/(27r) (dashed red) generated from foreground reduced WMAP full sky 
maps of W-band. 



each angular scale I to be distributed according to the Gaussian probability density 

, , 1 

Pi {aim) = '-'i (2) 

To evaluate the Gaussianity of an arbitrary distribution one usually looks at different statistical moments, which 
are quite noisy for a small number of data points. Since we are mainly interested in a statistics of the distribution 
tails |a;,„|^ > aCi it will be convenient to define^ 

a roc 

Cr = -— r^pi{r,6)2^rdr (3) 

1 + a Jv^i 

where we switched to polar coordinates {aim = re'^) and the normalization was chosen such that C" = Ci for a 
Gaussian (see Eq. On Fig. |[T|) we plot the smoothed Ci and generated from the foreground reduced WMAP 

full sky maps of W-band^ (the frequency which is known to have the lowest foreground contamination) . We choose a 



^ More generally, one could be interested to test Gaussianity in an arbitrary range aCi < |a;ml < f^Ci which is easily derived from C" 
and Cf . 

^ The foreground reduced maps were obtained from the Legacy Archive for Microwave Background Data Analysis website and the 
HEALPix C++ routines were used to generate the corresponding aj^'s. 
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Figure 2: Non-Gaussianity of foreground reduced WMAP full sky maps of W-band. 



value of a = 1 such that the presence of a non-Gaussianity peak is most pronounced and the results have a sufHcient 
statistical significance. One could clearly see an anomalous peak of the red dashed line in the range 35 ^ Z < 55. The 
feature can be better illustrated by binning and plotting only the deviations from Gaussianity. On Figs. ^ and ^ 
we plot Cf — Ci generated from two different frequency maps and see a very similar behavior in W- and V-bands. 
This indicates that unknown foregrounds are not likely to explain the presence of the non-Gaussianity peak. 

To understand the statistical significance of the result we first calculate an expected number of data points in the 
analyzed range at a given angular scale: 



N{1) = I 



Pi{r) 2'Kr dr — le 



Then the cosmic variance in the values of C" with N{1) data points^ is given by 



(5(0 



1 



Ci 



1 



(4) 



(5) 



which is approximately equal to 5{l) w 0-8-^ for a = 1. The small red error bars plotted on Figs. ^ and ^ represent 
this variance with binning taken into account. That would be all if we had a clean data set, but the foreground reduced 



We neglect possible variations in the number of data points which is justified for I > 20 and a < 1. 
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Figure 3: Non-Gaussianity of foreground reduced WMAP full sky maps of V-band. 

maps contain very high contamination form the galactic plane. The superposition of the galactic noise and CMB 
radiation increases significantly the total variance, which could be obtained from the observed variations at the angular 
scales away from the peak. The large blue error bars on Figs. ^ and ^ give a rough estimate of these uncertainties. 

To get a better insight, we reanalyze the data by penalizing the multipole moments which are most susceptible to 
the galactic noise. Clearly, the spherical harmonics 

Yira 0,6"'^'^ Pr{cos{0)), (6) 

with larger m (corresponding to larger relative amplitudes of Legendre functions near equator and more frequent 
oscillations of Yim along the azimuthal angle tp) pick up a stronger signal from the galaxy. To take into account 
unequal contamination of the multipoles we introduce weighting which increases the significance of relatively clean 
modes, 

w{aim) = m"'^, (7) 

where (3 > 0. On Fig. (JH) we plot non-Gaussianity of W-band map with weighting parameter f3 = I and error bars 
described by Eq. ([5]). One can see that the statistical variations associated with the galactic plane are greatly reduced, 
but the anomalous peak is still present. It follows that the probability for such an excessive non-Gaussianity to be 
present in a range from ^ « 36 to Z « 52 is about P « 0.000001.^ 



^ This result must be taken with caution, since the foreground reduced maps are highly contaminated at the galactic plane from where 
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Figure 4: Non-Gaussianity of CMB estimated by multipoles weighting of W-band map. 

The analysis of the WMAP ILC^ and Harmonic ILC''^ maps did not show any significant deviations from Gaussian- 
ity. Although, the linear combination techniques are known to introduce unphysical correlations between separate 
multipole moments, it remains to be seen whether this effect is significant enough to eliminate a non-Gaussian peak 
which is observed independently in different frequency bands. 

III. DISCUSSION 

A prehminary analysis of the WMAP foreground reduced maps revealed some unexpected non-Gaussianities around 
Z ^ 45 ± 10. The standard inflation does not predict such effects, but one might always attribute them to unknown 
astrophysical sources. It is quite likely that foreground contamination of the galactic plane would modify the tails of 
distributions, but it is not very likely that the same non-Gaussian feature is present at exactly the same location in 
different frequency bands* unless the feature has a cosmological origin. We would also like to note that the anomaly 
appears approximately where an outlying point of the WMAP Collaboration published power spectrum is located 



surprising non-Gaussianities migtit arise. 
^ It is well known that the WMAP ILC map is not reliable at small angular scales (/ > 20). 

We are grateful to Pavel Naselsky for providing us with the Harmonic ILC temperature map [l^ . 
* The Q-band foreground reduced map contains a very similar non-Gaussianity peak at the same angular scales, but this frequency is 

known to be most contaminated by foregrounds, and as a result the analysis of the Q-band map was not conclusive. 
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To check our results we analyzed further the W-band map, which contains the lowest contamination, but nevertheless 
quite noisy at the galactic plane. The main idea was to introduce unequal weighting of the multipole moments 
contaminated by foregrounds. The weighting technique allowed us to significantly reduce the uncertainties (See Fig. 
(0])) and to demonstrate that the observed non-Gaussianity is not likely to come from only galactic plane. Let us now 
describe two different mechanisms which could naturally lead to changing non-Gaussian tails at the angular scale of 
the anomaly. 

During inflationary stage the presence of non-scalar background fields (e.g. vector inflation could lead to a 
mixture of perturbations already at the Hnear order [H, Higher order corrections enhance the mixture (even 
after inflation leading to non-trivial correlations between Gaussian random variables (e.g. scalar and tensor 

components). Both effects combined could lead to appearance of non-Gaussian distribution tails. If we assume that 
primordial gravitational waves contribute a significant power to the CMB temperature fiuctuations and that the 
tensor modes became correlated with scalar modes, then changes in the behavior of one of the components could 
produce changes in non-Gaussianity. In fact the place where an anomalous peak appears is exactly where the tensor 
contribution to the temperature power spectrum should start a very rapid decline due to decaying evolution inside 
the horizon If this interpretation is correct, then one would expect the non-Gaussianity to change significantly 
close to the angular scales of ^ '^45. 

Although correlated gravitational waves could produce significant changes in non-Gaussianity, it is hard to argue 
that such scenarios are typical. In contrast, the production of cosmic strings is a generic feature of many models 
involving symmetry breaking phase transitions [l^. Non- linear interaction of the CMB photons with a network of 
cosmic strings could easily generate non-Gaussian fluctuations in the temperature. Even mild discontinuities of the 
CMB ternperature caused by cosmic strings are translated directly into appearances of non-Gaussian distribution 
tails [H, [131 ■ The angular scales of I ~ 200 correspond to the horizon size at the last scattering surface where the 
largest power from cosmic strings [HI, [131 and also from inflationary perturbation is expected. Thus, the relevant 
question is at which I the relative contribution from cosmic strings to the overall power is at its highest value. The 
answer to this highly non-trivial problem depends crucially on the evolution of the cosmic string network jH, [13, [13 
and in particular on the distribution of the cosmic string loops [13, [13, [13, [13, [131 ■ It is quite natural to expect that 
at Z ~ 45 where the power spectrum from inflationary perturbation just starts to rise, the relative contribution from 
cosmic strings is the highest. As a result, one should see the largest deviation from Gaussianity to be somewhat close 
to where the non-Gaussianity peak is in fact observed. 

Both mechanisms (cosmic strings or correlated gravitational waves) could in principle be employed to explain the 
anomaly, but before we seriously consider such scenarios it is important to fully understand the effect. This could 
only be achieved by performing a lot more sophisticated tests at the angular scales: 35 < ^ ^ 55. The detected non- 
Gaussianity could not have been seen from the usual bispectrum, but should in principle be observable through the 
trispectrum analysis (or through the forth statistical moment). It would also be desirable to probe larger multipoles 
{I > 100) in order to favor one model or the other. All of these issues deserve further investigation. 
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Note Added 

After the flrst version of the paper was submitted we learned from Pavel Naselsky another way to reduces the 
possible contamination from the galactic plane. Apparently the multipoles with odd I + m whose spherical harmonics 
vanish at the equator are less affected by the galactic noise. We had preformed a similar analysis with only odd I + m 
multipoles and discovered an opposite non-Gaussiniaty at the angular scales of the anomaly. Instead of having too 
many outlier we have to few, but the non-Gaussianity is still present at around Z 45 ± 10. It is hard to imagine 
how the galactic noise alone could act in opposite direction for odd and even I + m multipoles in only a narrow range 
of scales. In fact, similar non-Gaussianity was recently reported by another group [3l|. The authors looked at the 
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abundances of cold and hot spots in CMB and found a lack of power on the angular scales of several degrees. 
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